AD-A227  099 


NASA  Contractor  Report  181998 
ICASE  Report  No.  90-16 


ICASE 


SECONDARY  FREQUENCIES  IN  THE  WAKE  OF  A 
CIRCULAR  CYLINDER  WITH  VORTEX  SHEDDING 


Saul  S.  Abarbanel 
Wai  Sun  Don 
David  Gottlieb 
David  H.  Rudy 
Janies  C.  Townsend 


Contract  Nos.  NAS  1-18 107,  NAS  1-1 8605 
February  1990 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  Virginia  23665-5225 

Operated  by  the  Universities  Space  Research  Association 


NASA 

National  Aeronautics  and 
Space  Administration 


i 


Langley  Research  Center 

Hampton,  Virginia  23665-5225 


LiSTR i PUT; ON  STATEMENT'  A 

Anpro-'f-d  for  public  release; 
•V 'b'.’i-.o:'  Unlimited 


DTIC 


,  -a  ELECTE 
OCT  03  1990 


E 


Secondary  frequencies  in  the  wake  of  a  circular  cylinder  with 

vortex  shedding* 


Saul  S.  Abarbanel 
Tel  Aviv  University,  Tel  Aviv,  Israel 

Wai  Sun  Don 
David  Gottlieb 

Brown  University,  Providence,  RI,  USA 

David  H.  Rudy 
James  C.  Townsend 

NASA  Langley  Research  Center,  Hampton,  VA,  USA 


■  c<33ion  For 

LlIS  GPA&I 
.  ;  TAB 

L'v;-“-:ounced  □ 

•  Justification _ 


T1__ 

;  dv _ 

!  Distribution/ 

;  Availability  Codes 

1 

Avail  and/or 

;Di3t 

Spec 

lal 

Abstract 

A  detailed  numerical  study  of  two-dimensional  flow  past  a  circular  cylinder  at 
moderately  low  Reynolds  numbers  has  been  conducted  using  three  different  numerical 
algorithms  for  solving  the  time-dependent  compressible  Navier-Stokes  equations.  It 
was  found  that  if  the  algorithm  and  associated  boundary  conditions  were  consistent 
and  stable,  then  the  major  features  of  the  unsteady  wake  were  well-predicted.  However, 
it  was  also  found  that  even  stable  and  consistent  boundary  conditions  could  introduce 
additional  periodic  phenomena  reminiscent  of  the  type  seen  in  previous  wind-tunnel 
experiments.  However,  these  additional  frequencies  were  eliminated  by  formulating 
the  boundary  conditions  in  terms  of  the  characteristic  variables.  An  analysis  based  on 
a  simplified  model  provides  an  explanation  for  this  behavior. 
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Computer  Applications  in  Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton, 


VA  28665. 


1 


1.  Introduction 

Historically,  the  unsteady  wake  generated  by  a  circular  cylinder  in  low-speed  flow  has 
been  of  great  interest  to  computational  fluid  dynamicists  as  well  as  to  theoretical  and 
experimental  aerodynamicists.  The  Reynolds-number  range  between  40  and  1000  has 
been  of  particular  interest  because  it  spans  the  transition  from  steady  flow  to  unsteady 
wake  flow  dominated  by  the  periodic  shedding  of  vortices  from  the  cylinder.  The  shedding 
frequency  of  these  vortices  increases  with  Reynolds  number  over  this  range,  asymptotically 
approaching  a  constant  value.  However,  Sreenivasan  (1985)  measured  more  than  one 
distinct  frequency  in  the  shedding  regime  at  low  Reynolds  numbers.  In  addition  to  the 
vortex-shedding  frequency,  he  found  clearly  discernible,  lower  frequencies  in  the  frequency 
spectrum  for  the  streamwise  velocity  measured  in  the  wake.  These  additional  frequencies 
were  not  subharmonics  of  the  primary  shedding  frequency.  He  concluded  that  this  was  a 
feature  of  the  initial  stages  of  transition  to  turbulence  in  agreement  with  the  route  to  chaos 
described  by  Ruelle  and  Takens  (1971).  Sirovich  (1985)  suggested  that  these  additional 
modes  of  oscillation  could  be  described  theoretically  in  terms  of  the  classical  von  Karman 
vortex  street.  On  the  other  hand,  based  on  measurements  of  vibrating  and  nonvibrating 
cylinders,  van  Atta  and  Gharib  (1987)  concluded  that  the  additional  frequencies  found  by 
Sreenivasan  were  due  to  the  aeroelastic  coupling  of  the  vortex  wake  with  cylinder  vibration 
modes.  For  a  nonvibrating  cylinder,  they  found  no  spectral  peaks  other  than  the  primary 
Strouhal  shedding  frequency.  Recently,  however,  Sreenivasan  (1989)  determined  that  the 
cylinders  did  not  vibrate  in  his  previous  experiments. 

Conflicting  results  have  also  come  from  the  numerical  studies  of  this  flow.  Karniadakis 
and  Triantafyllou  (1989)  found  no  secondary  frequency  in  their  computation  of  the  incom¬ 
pressible  flow  past  a  circular  cylinder  using  the  spectral-element  method.  Moreover,  they 
were  able  to  excite  a  secondary  mode  by  introducing  an  external  forcing  function  into  the 
momentum  equations.  Townsend  et  al.  (1987)  did  find  a  secondary  frequency  in  their 
finite-difference  computation  of  low-speed  compressible  flow  past  a  circular  cylinder.  This 
low  frequency  was  found  in  the  frequency  spectrum  of  the  pressure  at  various  points  in 
the  wake.  However,  the  frequency  was  dependent  upon  the  size  of  the  solution  domain, 
suggesting  that  it  might  have  been  a  numerical  effect. 

The  present  paper  describes  the  results  of  further  numerical  studies  which  demonstrate 
that  the  secondary  frequency  found  in  the  compressible-flow  calculations  was  produced  by 
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the  far-field  boundary  conditions.  Three  different  numerical  methods  were  used  to  solve 
the  time-dependent  compressible  Navier-Stokes  equations.  The  first  of  these  was  the  finite- 
difference  technique  used  by  Townsend  et  al.  ( 1987) ,  the  second  method  was  a  fully-spectral 
method  (Don  1989),  while  the  third  was  a  mixed  spectral-finite-difference  method  (Don 
1989)  which  was  a  combination  of  the  other  two  methods.  Computations  were  made  with 
the  finite-difference  code  with  two  sets  of  far-field  boundary  conditions.  These  calculations, 
which  were  made  for  flow  conditions  in  both  the  vortex-shedding  regime  and  the  steady- 
wake  regime,  demonstrate  the  effect  of  the  boundary  conditions  on  the  frequency  spectrum 
of  the  pressure  in  the  wake.  With  the  proper  choice  of  far-field  boundary  conditions, 
no  frequencies  other  than  the  primary  shedding  frequency  and  multiples  of  the  primary 
frequency  were  found  at  a  Reynolds  number  of  80.  In  addition,  when  the  “improper” 
boundary  conditions  were  used  to  compute  the  flow  at  a  Reynolds  number  of  20,  where  the 
flow  is  known  to  be  steady,  then  indeed  no  shedding  frequency  was  predicted  although  the 
secondary  frequency  did  appear  in  the  solution.  These  findings  were  further  substantiated 
by  calculations  with  the  highly-accurate  spectral  method. 

The  details  of  the  solution  procedures,  including  boundary  conditions,  are  described  in 
§2  of  the  paper.  All  of  the  computed  results  are  presented  in  §3.  In  §4,  a  mathematical 
analysis  of  boundary  conditions  for  a  model  system  is  presented  which  shows  how  additional 
frequencies  can  be  produced  by  the  boundary  conditions.  Finally,  a  summary  of  the  study 
and  some  concluding  remarks  are  given  in  §5. 

2.  Problem  formulation  and  numerical  methods 

2.1  Governing  equations 

For  the  present  compressible- flow  calculations,  the  governing  equations  are  those  which 
describe  the  conservation  of  mass,  momentum,  and  energy  of  an  ideal  fluid  in  the  absence 
of  external  forces.  The  nondimensionalform  of  these  equations  in  a  general  curvilinear 
coordinate  system  is 
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J  is  the  Jacobian  of  the  transformation  from  Cartesian  coordinates  (x,y)  to  the  general 
curvilinear  coordinates  (£,77), 
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and  the  specific  internal  energy,  E,  is  defined  by 


E  —  p 

The  elements  of  the  stress  tensor  are 
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The  nondimensional  equation  of  state  is 


(2.6) 


p  =  h  - 1 )pT, 


(2.7) 


and  the  viscosity  was  given  by  the  Sutherland  law 


H  = 


C{T  § 
C2  +  T 


(2.8) 


where  C\  and  C2  are  constants.  The  ratio  of  specific  heats  7  was  1.4,  and  the  Prandtl  num¬ 
ber  Pr  was  0.72.  These  equations  were  nondimensionalized  using  the  following  reference 
quantities:  Uref  =  U^,  pref  =  /><»,  pTt /  =  PocU^,  and  Tref  =  H£/cu  where  and  p^  are 
the  free-stream  velocity  and  density,  respectively,  and  cv  is  the  specific  heat  at  constant 
volume.  Therefore,  the  reference  Reynolds  number  is  defined  as  Reref  ■=  UrefD/vref  where 
D  is  the  diameter  of  the  cylinder  and  i/ref  is  the  kinematic  viscosity  based  on  the  reference 
temperature.  However,  it  should  be  noted  that  the  Reynolds  number  used  in  subsequent 
sections  of  the  paper  to  characterize  the  flow  conditions  for  which  calculations  were  made 
is  based  upon  free-stream  quantities,  i.e.,  —  U^D/u^. 


2.2  Grid 


A  polar  grid  was  used  in  the  calculations.  The  grid  was  generated  in  the  physical  (x,y) 
plane  and  mapped  numerically  onto  the  computational  (£,7)  plane.  The  £  coordinate 
corresponded  to  the  circumferential  direction  and  the  r)  coordinate  to  the  radial  direction. 
The  grid  was  stretched  in  both  directions  so  that  grid  points  were  clustered  toward  the 
cylinder  in  the  radial  direction  and  toward  the  wake  region  in  the  circumferential  direction. 
For  the  finite-difference  calculations,  the  basic  grid  had  122  points  in  the  circumferential 
direction  and  151  points  normal  to  the  body.  The  outer  boundary  of  the  grid  was  located  25 
diameters  from  the  cylinder  surface.  For  the  computations  with  the  spectral  method,  the 
grid  had  64  and  48  points  in  the  angular  and  radial  directions,  respectively.  Calculations 
were  made  with  the  outer  boundary  at  both  20  and  22.5  diameters  away  from  the  cylinder. 
For  computations  with  the  mixed  spectral-finite-difierence  method,  the  grid  had  70  and 
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100  points  in  the  circumferential  and  radial  directions,  respectively.  For  this  case,  the  far- 
field  boundary  was  located  22.5  diameters  from  the  cylinder  surface.  In  the  experiments 
of  Sreenivasan  (1985),  the  wind-tunnel  walls  were  located  30  diameters  above  and  below 
the  cylinder1. 


2.3  Numerical  methods 

The  finite-difference  method  used  in  the  present  study  was  the  original  unsplit  tech¬ 
nique  of  MacCormack  (1969),  which  has  second-order  accuracy  in  both  space  and  time. 
The  method  is  an  explicit,  conditionally-stable,  predictor-corrector  scheme.  Forward  dif¬ 
ferences  were  used  to  approximate  the  derivatives  of  the  fluxes  in  equation  (2.1)  in  the 
predictor  step,  and  backward  differences  were  used  in  the  corrector  step.  The  first  deriva¬ 
tives  appearing  in  the  viscous  flux  terms,  Fv  and  Gv,  were  approximated  with  backward 
differences  if  the  flux  derivative  was  being  approximated  with  a  forward  difference,  and 
vice-versa. 

The  spectral  algorithm  used  both  the  Fourier  and  Chebyshev  collocation  methods. 
Because  of  the  polar  grid,  the  flow  could  be  treated  as  being  periodic  at  the  boundaries 
of  the  solution  domain  in  the  circumferential  direction.  Thus,  the  Fourier  collocation 
(pseudospectral)  method  was  the  natural  choice  for  the  circumferential  direction.  Since 
the  flow  in  the  radial  direction  was  not  periodic,  the  Chebyshev  collocation  method  was 
used  in  that  direction.  Both  of  these  methods  were  employed  in  the  form  of  the  matrix 
vector  multiplication  method  instead  of  the  more  commonly  used  Fast  Fourier  Transform 
method.  To  improve  the  stability  of  the  algorithm,  a  fourth-order  exponential  filter  (Don 
1989)  was  applied  to  the  differentiation  and  solution  matrices.  The  solution  was  advanced 
in  time  using  a  second-order  Runge-Kutta  method. 

A  third  approach  which  combined  these  two  methods  was  also  used.  In  this  case,  the 
Fourier  method  was  used  in  the  circumferential  direction  and  the  finite-difference  method 
in  the  radial  direction. 


2.4  Boundary  conditions 

A  sketch  of  the  solution  domain  in  the  physical  plane  is  shown  in  figure  1.  Periodic 
boundary  conditions  for  the  dependent  variables  were  applied  along  the  cut  line  in  the 
‘Private  communication. 
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wake  for  both  the  finite-difference  and  spectral  methods.  Standard  solid-wall  boundary 
conditions  for  viscous  flow  were  applied  at  the  cylinder  surface.  A  no-slip  condition  was 
applied  so  that  the  velocity  components,  u  and  v,  were  specified  to  be  zero  at  the  surface. 
The  wall  temperature  was  held  fixed  at  the  value  of  the  free-stream  temperature.  In 
the  finite-difference  code,  the  density  at  the  wall  was  obtained  from  an  extrapolation  of 
interior-point  values  of  pressure  in  the  direction  normal  to  the  wall.  In  the  spectral  code, 
the  density  at  the  wall  was  computed  as  a  part  of  the  solution. 

The  outer  boundary  is  an  arbitrarily-chosen  boundary,  introduced  only  to  restrict  the 
computational  domain  to  a  finite  size.  Since  this  is  a  boundary  across  which  fluid  passes 
either  into  or  out  of  the  computational  region  and  since  disturbances  can  propagate  up¬ 
stream  as  well  as  downstream  in  a  subsonic  flow,  careful  consideration  must  be  given  to 
the  specification  of  boundary  conditions  on  this  outer  boundary.  If  this  outer  boundary  is 
located  sufficiently  far  from  the  cylinder  surface,  viscous  effects  axe  negligible  in  the  flow 
crossing  the  boundaries  except  in  the  narrow  region  where  the  cylinder  wake  is  located. 
As  a  result,  the  proper  choice  of  boundary  conditions  can  be  found  from  an  analysis  of  the 
inviscid  form  of  the  governing  equations. 

For  a  subsonic  inflow  boundary,  this  analysis  shows  that  there  are  three  characteristics 
coming  into  the  solution  domain  and  one  outgoing  characteristic.  Thus,  there  must  be 
three  quantities  specified  at  an  inflow  boundary.  Since  there  are  four  governing  equations, 
the  numerical  method  requires  a  fourth  boundary  condition  in  addition  to  the  three  re¬ 
quired  for  proper  specification  of  the  boundary  conditions  in  a  mathematical  sense.  In  the 
finite-difference  code,  calculations  were  made  with  two  different  sets  of  inflow  boundary 
conditions.  The  first  set,  which  was  used  in  the  previous  solutions  reported  by  Townsend 
et  al.  (1987),  was  formulated  in  terms  of  the  primitive  variables  u,v,p,  and  T.  At  the 
inflow  boundary,  free-stream  values  of  the  two  velocity  components  and  the  temperature 
were  specified,  i.e., 


W inflow  —  Moo 

(2.9) 

V inflow  ^oo 

(2.10) 

Tin  flow  =  Too 

(2.11) 

The  density  was  obtained  from  a  zeroth-order  extrapolation  of  pressure  from  the  interior 
of  the  solution  domain  and  the  use  of  equation  (2.7).  In  preliminary  calculations  in  whirK 
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only  the  primary  shedding  frequency  was  of  interest,  this  boundary  condition  gave  results 
which  were  not  significantly  different  from  those  obtained  using  an  extrapolation  of  the 
outgoing  characteristic  to  obtain  pressure. 


At  the  outflow  boundary,  only  one  analytic  boundary  condition  is  needed  since  there  is 
only  one  incoming  characteristic.  A  typical  choice  is  the  specification  of  the  static  pressure. 
However,  as  demonstrated  numerically  by  Rudy  and  Strikwerda  (1980),  this  can  cause 
waves  to  be  reflected  from  the  outflow  boundary  back  into  the  solution  domain,  adversely 
affecting  the  solution  in  the  interior  of  the  domain.  A  nonreflecting  boundary  condition 
based  on  the  work  of  Engquist  and  Majda  (1977)  was  used  in  the  present  study.  Thus,  the 
pressure  at  the  outflow  boundary  was  found  from  a  finite-difference  approximation  to  the 
equation 


dp 


(  du 


dv 


-57  -Pc\  u°o  "XT'  =  0 


dt 


^  dt 


dy, 


where  c  is  the  nondimensional  speed  of  sound  given  by 


(2.12) 


c  =  -  1  )T  , 


(2.13) 


This  boundary  condition  was  applied  along  the  outflow  boundary  where  the  viscous  wake 
crossed  the  boundary  as  well  as  the  region  immediately  above  and  below  the  wake.  Along 
the  remainder  of  the  outflow  boundary,  where  the  flow  was  essentially  inviscid,  the  pressure 
was  specified  to  be  the  free-stream  pressure.  The  variables  u,v,  and  T  along  the  entire 
outflow  boundary  were  obtained  from  zeroth-order  extrapolation,  and  the  density  was  then 
obtained  from  equation  (2.7)  using  the  boundary  value  of  p  and  the  extrapolated  value  of 


T. 


The  second  set  of  boundary  conditions,  which  was  also  used  in  the  spectral  calculations, 
was  based  entirely  on  the  characteristic  variables.  An  analysis  of  the  inviscid  form  of 
equation  (2.1),  linearized  around  the  free-stream  conditions,  shows  (Gottlieb,  Lustman, 
and  Streett  1984)  that  the  characteristic  variables,  with  corresponding  eigenvalues  ax  = 
a2  =  U  •  N ,  a3  =  U  •  N  —  c,  and  a4  =  U  •  N  +  c,  are 

Ri  =  -  (jpUoo  ■  tfoo  -  M  •  C/oo  +  tf)  (2.14) 

Ri  =  (mu  -  pu^rjy  -  (m„  -  pv^ijx  (2.15) 

Rz  =  -(Af  -  pUoo)  •  N  +  (~~)  •  Uoo  ~  M  ■  Uoo  +  (2.16) 

R<  =  (Af  -  pt/oo)  •  N  +  (ip!)  Qpt/oo  •  -  M  •  f/oo  +  #)  (2.17) 
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where  —  (u^,^)  is  the  free-stream  velocity,  M  =  (mu,mu)  =  (pu,pv)  is  the  local  mo¬ 
mentum,  Coo  is  the  free-stream  speed  of  sound,  and  N  =  ( rjx ,  rjy )  =  (|^,  |^)/ ^(fl)2  +  (f^)2 
is  the  unit  outward  normal  vector.  At  the  subsonic  inflow  boundary,  the  characteristic  vari¬ 
ables  corresponding  to  the  three  incoming  characteristics  were  specified  using  free-stream 
values.  These  boundary  values  then  became 

F,  =  ^(PocMoo^oo),  t  =  l,2,3  (2.18) 

Furthermore,  the  required  numerical  boundary  condition  is 

F4  —  Rii'Pnumt  A/fium>  Fntim)  (2.19) 

where  values  in  F4  are  obtained  by  extrapolation  from  the  interior  of  the  domain  in  the 
finite-difference  code  and  are  obtained  as  a  part  of  the  solution  in  the  spectral  method. 
Similarly,  at  the  outflow  boundary,  the  characteristic  variable  corresponding  to  the  incom¬ 
ing  characteristic  was  specified,  i.e., 

Fz  =  Rz{poo,M00,E!X>)  (2.20) 

and  the  values  for  the  other  three  characteristic  variables  were  found  numerically  as  at  the 


inflow  boundary,  so  that 

Fi=  Ri{pnum,Mnum,Enum),  i  =  1,2,4.  (2.21) 

At  each  boundary,  the  conserved  variables  were  then  computed  from  the  appropriate  sys¬ 
tem  of  equations,  giving 

p  =  — (Fj  +  4>i)  (2.22) 

Loo 

mu  =  "f”  ^3^7y)  "f“  Puoo  (2.23) 

mv  =  {<)> 2Vv  -  <f> 3*7z)  +  pv 00  (2-24) 

E  ~  (j)  1  -—pUoQ  •  U 00  +  muUoo  +  m^oo  (2.25) 

& 

on  the  outer  boundary,  where 


«  _  ^oo  {F3  "h  F4)  _  ip\  _  v 

4>i  -  - - - - - ,  <t>2  —  x(F4  -  Fa  ,  4>3  —  F2  . 

7  -  1  2  2 
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3.  Results  and  discussion 

Computations  were  made  with  all  three  methods  for  Mach  0.4  flow  at  a  Reynolds 
number  of  80.  The  initial  flow  field  was  specified  to  be  free-stream  flow  everywhere  except 
at  the  cylinder  surface.  Figure  2(a)  shows  pressure  contours  (isobars)  from  the  finite- 
difference  solution  using  the  primitive-variable  far-field  boundary  conditions  at  a  point  in 
time  after  the  periodic  vortex  shedding  had  been  fully  established.  The  pressure  field  for 
finite-difference  solution  using  the  characteristic  far-field  boundary  conditions  is  shown  in 
figure  2(b).  The  corresponding  solution  for  the  fully-spectral  code  is  given  in  figure  3. 

As  the  computations  were  made,  the  values  of  the  pressure  were  saved  at  ten  selected 
grid  points  in  the  flow  field.  At  least  four  of  these  grid  points  were  located  in  the  wake 
region  and  the  others  were  upstream  of  the  cylinder.  For  the  finite-difference  code,  the 
computed  pressure  data  were  saved  every  twenty  time  steps,  and  for  the  fully-spectral 
and  mixed  spectral-finite-difference  codes,  because  of  the  larger  time  step  used  in  these 
methods,  the  data  were  saved  after  every  time  step. 

The  set  of  pressure  data  at  each  point  was  analyzed  for  its  frequency  spectrum  using 
a  discrete  Fast  Fourier  Transform  (FFT)  technique.  First,  the  time-averaged  pressure  was 
subtracted  from  the  pressure  data  to  form  a  time  sequence  pj  of  N  points,  where  N  =  2m 
for  some  positive  integer  m.  Second,  the  Fourier  coefficients  Fk  of  the  transformation 

^v-"1  j(  2rjk  i 

Fk  =  A  ^2  Pie  '  N  '  for  k  =  0, 1, . . . ,  N  —  1  (3.1) 

j=o 

were  obtained  using  a  pre-coded  FFT  subroutine.  A  =  1  in  the  subroutine  used  to  analyze 
the  finite-difference  results  and  A  =  l/N  in  the  subroutine  used  to  analyze  the  results 
from  the  fully-spectral  and  mixed-method  codes.  The  power  spectral  density  dk  for  wave 
number  k  was  found  then  simply  as 

dk  =  B\Fk\2  where  k  =  0, 1, . . .  ,jV  —  1.  (3.2) 

B  —  2/N  for  the  finite-difference  results  and  B  =  1  for  the  others.  The  Strouhal  number 
(nondimensional  frequency)  corresponding  to  the  vortex-shedding  frequency  was 

St1  =  ki/{NAt),  (3.3) 

where  At  was  the  time  step  and  kx  was  the  value  of  k  such  that 

dki  =  max  dk  . 
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Similarly,  St2  is  the  Strouhal  number  corresponding  to  the  secondary  frequency. 


Note  that,  because  the  power  density  was  found  only  for  discrete  frequencies,  there 
was  an  inherent  discretization  error  in  determining  the  dominant  frequencies  in  the  flow. 
The  actual  peak  frequency  would  be  within  half  of  the  discretization,  or  ±l/(2iVA<),  of 
that  given  by  the  FFT.  In  some  of  the  present  cases,  the  resulting  error  in  the  Strouhal 
number  could  be  as  large  as  0.003.  In  an  effort  to  get  a  somewhat  better  estimate,  for  each 
dominant  wave  number  a  peak  in  the  frequency  spectrum  was  estimated  using  a  three- 
point  parabolic  interpolation.  These  values  were  used  in  table  1,  which  summarizes  the 
results  found  in  the  present  study. 


Method 

tD 

Moo 

Domain 

Inflow  BC 

Outflow  BC 

St  i 

st2 

FD 

mm 

51D 

Primitive 

Primitive 

FD 

m 

51D 

Primitive 

Primitive 

1 

.0347 

FD 

80 

51D 

Primitive 

Primitive 

EM 

FD 

80 

0.4 

51D 

Characteristic 

Characteristic 

.1518 

— 

S 

80 

0.4 

46D 

Characteristic 

Characteristic 

.1574 

— 

S 

80 

0.4 

41D 

Characteristic 

Characteristic 

.1565 

— 

M 

80 

0.4 

46D 

Primitive 

Primitive 

.1589 

.0225 

M 

80 

idj 

46D 

Characteristic 

Characteristic 

.1558 

— 

— 

Table  I:  Shedding  frequencies  computed  with  finite-difference  code  (FD),  mixed  fi¬ 
nite-difference-spectral  code  (M),  and  fully-spectral  code  (S). 


These  frequencies  were  based  on  the  computed  pressure  at  a  point  in  the  wake  10  di¬ 
ameters  downstream  of  the  cylinder  and  one  diameter  above  the  wake  centerline.  This 
location  corresponds  to  the  principle  one  used  in  the  experiments  of  Sreenivasan  (1985).  It 
should  be  noted  that  the  frequencies  for  the  finite-difference  solution  with  primitive  vari¬ 
ables  at  Mach  0.4  differ  slightly  from  those  given  by  Townsend  et  al.  (1987),  in  which  the 
results  were  incorrect  due  to  a  misinterpretation  of  the  output  firm  the  FFT  subroui’ne. 

Experimental  measurements  of  the  primary  shedding  frequency  in  the  wake  of  a  circular 
cylinder  at  low  Reynolds  numbers  were  reported  by  Roshko  (1953).  An  approximation  to 
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I 

the  data  was  given  by  the  equation 

(21  2  \ 

•  (3-4) 

-flCoo tD  ) 

For  a  Reynolds  number  of  80,  this  equation  yields  S£i=0.1558.  As  shown  in  table  1,  the 
computed  values  of  St i  are  very  close  to  the  experimental  value  in  all  cases.  The  largest 
deviation  is  3.47%  and  most  of  the  values  are  within  2%  of  the  experimental  frequency. 
Therefore,  the  primary  vortex-shedding  frequency  in  the  wake  is  well-predicted  with  either 
set  of  boundary  conditions. 

Plots  of  the  time  history  of  the  pressure  and  the  resulting  frequency  spectrum  at  the 
selected  point  in  the  wake  flow  from  the  finite-difference  calculations  are  shown  in  figure  4 
for  both  sets  of  boundary  conditions.  The  corresponding  plots  for  the  fully-spectral  calcu¬ 
lations  with  characteristic  boundary  conditions  are  shown  in  figure  5.  As  shown  in  figure  4, 
the  solution  with  primitive-variable  boundary  conditions  produces  a  low  frequency  which 
is  not  seen  in  the  spectrum  from  the  calculations  with  the  characteristic  boundary  condi¬ 
tions.  This  frequency  is  not  a  subharmonic  of  the  primary  shedding  frequency  and  was 
not  found  in  any  of  the  solutions  using  any  of  the  three  codes  when  characteristic  bound¬ 
ary  conditions  were  employed.  Furthermore,  the  amplitude  of  the  secondary  frequency 
remained  relatively  constant  at  all  locations  for  which  the  spectrum  was  computed,  in¬ 
cluding  a  location  10  diameters  upstream  of  the  cylinder.  At  this  point,  the  amplitude 
of  the  primary  shedding  frequency  was  significantly  lower  than  in  the  wake  so  that  the 
secondary  frequency  was  the  dominant  frequency  in  the  signal.  These  factors  suggest  that 
the  secondary  frequency  is  a  spurious  numerical  artifact. 

Calculations  were  also  made  with  the  finite-difference  code  using  the  primitive-variable 
boundary  conditions  to  determine  the  effect  of  Mach  number  on  the  computed  frequencies. 
Computations  were  made  at  Mach  numbers  of  0.2  and  0.1  for  a  Reynolds  number  of  80. 
As  shown  in  table  1,  there  was  only  a  small  effect  on  St i  as  Mach  number  was  lowered. 
However,  as  Mach  number  was  decreased  by  a  factor  of  2,  St-z  increased  by  a  factor  of 
approximately  2. 

To  further  establish  that  the  secondary  frequency  is  a  numerical  artifact,  calculations 
were  made  for  Mach  0.4  flow  at  a  Reynolds  number  of  20.  Under  these  conditions,  the  flow 
in  the  wake  should  be  steady,  consisting  of  two  symmetric  counterrotating  vortices  just 
behind  the  cylinder.  Calculations  with  the  finite-difference  code  produced  the  expected 
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steady  wake  with  both  sets  of  boundary  conditions.  However,  when  the  primitive-variable 
boundary  conditions  were  used,  a  secondary  frequency  appeared  in  the  frequency  spectrum 
of  the  pressure  in  the  wake  with  a  value  of  5^2=0.0158.  This  spurious  secondary  frequency 
did  not  appear  when  characteristic  boundary  conditions  were  used. 

In  summary,  the  calculations  presented  in  this  section  demonstrate  that  a  secondary 
frequency  can  be  introduced  into  the  solution  by  the  far-held  boundary  conditions.  By 
coincidence,  the  value  of  this  frequency  was  very  close  to  that  found  experimentally  by 
Sreenivasan  (1985).  This  frequency  was  not  found  when  the  boundary  conditions  are 
properly  formulated  in  terms  of  characteristic  variables. 

4.  Analysis 

As  shown  in  §3,  the  secondary  frequency  disappeared  when  boundary  conditions  based 
on  the  characteristic  variables  were  used  at  the  outer  boundary,  demonstrating  that  the 
secondary  frequency  computed  with  the  primitive-variable  boundary  conditions  was,  in 
fact,  of  numerical  origin.  In  this  section,  a  theoretical  explanation  is  given  for  these  results 
using  a  modal  analysis  of  the  effect  of  inflow  boundary  conditions  on  the  stability  of  small 
perturbations  to  the  flow.  It  is  shown  that  the  use  of  boundary  conditions  based  on 
noncharacteristic  variables  causes  a  temporally-periodic  perturbation.  Furthermore,  it  is 
demonstrated  that  the  use  of  characteristic  variables  eliminates  this  periodic  disturbance. 

The  theoretical  study  is  based  upon  a  one-dimensional  treatment  of  the  Euler  equations. 
The  use  of  the  inviscid  form  of  the  governing  equations  is  justified  at  distances  far  from 
the  cylinder  since  the  Reynolds  numbers  under  consideration  place  the  flow  well  outside 
the  limits  of  the  Stokes  and  Oseen  regions.  The  assumption  of  one-dimensionality  is  more 
problematical.  However,  it  is  a  reasonably  good  representation  for  the  calculations  in  which 
a  finite-difference  method  was  used  in  the  radial  direction  and  the  Fourier  collocation 
method  was  used  in  the  circumferential  direction.  In  fact,  it  will  be  shown  that  the 
theoretical  predictions  of  the  secondary  frequency  are  in  good  agreement  (within  15%) 
with  the  computed  values  for  this  case.  The  agreement  with  the  computed  results  from 
the  fully  finite-difference  code  is  not  as  good,  although  the  trend  with  Mach  number  is 
well-predicted. 

The  analysis  considers  a  region  of  flow  along  the  x  axis  sufficiently  far  upstream  from 
the  cylinder.  The  one-dimensional  form  of  equation  (2.1),  under  the  assumption  of  inviscid 


12 


flow  then  becomes 


dU  dF  n 
aT  +  ^-°' 


(4.1) 


where 


U  = 


P 

m 

E 


F  = 


m 

T  +  P 


[  U(E  +  p)  j 

and  m  —  pu.  In  terms  of  the  dependent  variables  in  U,  the  equation  of  state,  equation 
(2.7),  is 


(7-1) 


l.m\ 

E  -  -  — ) 

2k  p  ’ 


Equation  (4.1)  can  be  rewritten  in  nonconservation  form  as 

dU  *,mdU  n 

-  +  A(J/)-=0, 


(4.2) 


(4.3) 


_  dF 
dU 

Linearizing  about  the  steady  free-stream  conditions,  this  equation  becomes 


w’ere  A(U)  =  |^  is  the  Jacobian  of  the  flux  with  respect  to  the  solution  vector  U. 


+  A(U~)A(6U)  =  0, 


(4.4) 


where  6U  =  U  —  Uoo  is  the  perturbation  vector.  The  matrix  A(Uoo)  =  ■  ^  '  _  has 

J  U  —  Uoo 

three  eigenvalues  d\  —  Uoo  —  Coo,  a 2  =  «oo  +  Coot  and  03  =  u The  corresponding 
eigenvectors  in  terms  of  the  conserved-variable  perturbations  are 


Ri  =  -(6m  -  Uoo 6p)  + 

R2  =  (6m  -  Uoo6p)  +  - 


/y  _  11 

-  -u^p  -  u006m  +  6E 

Coo  La 

\ulob  ~  Uoo6m  +  6E 


and 


—  Cqq6  p 


7-1  ri 


1  2  c 

2Uoo 6p  -  «0 


,6m  +  6E 


(4.5) 

(4.6) 

(4.7) 
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Furthermore,  using  a  linearized  version  of  the  equation  of  state,  i.e., 

1 


Bp  =  (7  -  1) 


l2U°°Sp  ~  u^m  +  BE 


(4.8) 


equations  (4.5)-(4.7)  may  also  be  written  in  terms  of  the  primitive-variable  perturbations 
(Bp,Bu,Bp),  giving, 


R\  =  -PooSu  H - Bp, 

Coo 

Rl  =  PooBu  -I - Bp, 

Coo 


and 


Rl  —  CooBp 


—Bp. 

Coo 


(4.9) 

(4.10) 

(4.11) 


Any  perturbation  imposed  on  the  free-stream  solution  will  evolve  as  a  combination  of  these 
eigenvectors.  This  evolution  can  be  studied  locally  near  the  outer  computational  boundary 
by  employing  the  modal  analysis  developed  by  Gustafsson,  Kreiss,  and  Sundstrom  (1972) 
and  Osher  (1969)  to  study  the  stability  of  numerical  approximations  to  initial-boundary- 
value  problems. 

Consider  an  inflow  boundary  point  with  subsonic  flow  crossing  the  boundary  into  the 
solution  domain.  At  this  point,  <12  =  «oo  +  >  0  and  a3  =  «<»  >  0-  Therefore,  the 

characteristic  variables,  R2  and  R$,  corresponding  to  these  two  incoming  characteristics 
must  be  specified.  The  third  one,  Ri,  corresponds  to  ax  =  <  0  and  thus  exits 

the  domain.  To  model  the  situation  near  the  inflow  boundary,  consider  the  governing 
equations  in  terms  of  the  characteristic  variables,  viz., 


or 


dRi  .  ,  dRi 

at  +  (“°°  ax 

=  0, 

(4.12) 

dRi  ,  ,  ,  ^dR2 

dt  +  («o, 0  +  Coo)  dx 

=  0, 

(4.13) 

dR3  dRi 

(4.14) 

dt  +  Uo°  dx 

=  0, 

dR,  dR. 

(4.15) 

at  +a‘  ax  =  °-  s 

=  1,2,3. 

In  the  one-dimensional  case,  the  finite-difference  algorithm  described  in  §3  is  equivalent 
to  a  Lax-Wendroff  scheme  and  can  be  written  as 

»,n+l  _  OjAt  f „„s,n  \  ,  ®»(Af) 


W 


9>n  0>8LXl  /  a  n  a,n  \  .  aal^V  (  a,n  0  a,n  ,  a,n  \  •  .  «  i  A  ^ 

=  j  ~  TKi  V^Vi  ~  *i-x)  +  K+i  -  2wi  +  wi-i) >  J  >  °>  (4-16) 
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where  w‘,n  =  w‘(j  Ax,nAt)  is  the  finite-difference  approximation  to  Rs(x,t).  Equation 
(4.16)  is  solved  with  inflow-type  boundary  conditions  imposed  at  Xb0Undary  =  (iAi)J=0 
which  are  of  the  form 

(i?2  —  aRi)j=0  =  0, 

(i?3  -  0Ri)j=o  —  0, 

and 

(i?i  +  0R2  +  eRs)j=0  —  [Ri  +  0R2  +  ei?3)J=1.  (4.19) 

Thus,  equation  (4.19)  is  a  zeroth-order  extrapolation  on  Rx. 

The  modal  solution  ansatz  for  this  case  is 

Ra  =  A„K\zn  s  =  1,2,3.  (4.20) 

Substituting  this  solution  into  equation  (4.16)  yields  the  three  modal  characteristic  equa¬ 
tions  given  by 

(A?  -  As)/cs2  +  2(1  -  ^  -  A?)>c,  +  (A*  +  A,)  =  0,  a=  1,2,3  (4.21) 

where  A,  =  2^.  Substituting  equation  (4.20)  into  equations  (4.17)  to  (4.19)  gives  the 
modal  representation  of  the  boundary  conditions, 

A2  =  aAlt  {Ai  ±  0)  (4.22) 

A3  =  (3A\  (4.23) 

and 

Ai(l  -  fCi)  +  crA2(l  -  k2)  +  tAz(  1  -  /c3)  =  0.  (4.24) 

Since  A\  ^  0,  equation  (4.24)  may  be  rewritten  as 

(1  -  Ki)  +  ao(l  -  /c2)  +  U  ~~  *3)  =  0.  (4.25) 

It  is  necessary  to  find  z  such  that  |  z  |=  1  and  k„  (s  =  1,2,3)  such  that  |  k,  |<  1  which 
solve  equations  (4.21)  and  (4.25)  simultaneously.  It  should  be  noted  that  |  z  |=  1,  z  ^  l, 
corresponds  to  a  solution  that  is  purely  periodic  in  time.  If  the  search  for  such  solutions  fails 


(4.17) 

(4.18) 
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(which  it  does,  as  it  will  be  shown,  when  a  —  (3  =0,  i.e.,  for  characteristic  specification), 
then  the  boundary  treatment  does  not  introduce  any  spurious  frequency.  If,  under  a  given 
set  of  boundary  conditions,  there  is  a  solution  to  equations  (4.21)  to  (4.25),  then  the 
phase  of  a  gives  the  temporal  frequency  of  the  numerically-introduced  perturbation.  Of 
course,  it  is  hoped  that  at  least  one  of  the  k’s  will  be  nearly  1  in  magnitude.  Otherwise, 
the  temporally-periodic  perturbation  will  decay  spatially  as  j/cj-? .  In  fact,  in  all  of  the 
calculations  with  the  model,  it  is  al  vays  the  case  that  .995  <j  K2  J<  1. 

In  the  actual  numerical  computations  presented  in  §3,  the  noncharacteristic  boundary 
conditions  (translated  to  the  present  one-dimensional  model)  consisted  of  specifying  the 
velocity  and  temperature  ((<5u)y-o  —  0,  (/>T);_0  =  0)  and  performing  the  extrapolation 

=  {p)}=i^~-  (4.26) 

CO 

This  is  equivalent  to 

(H=o  =  (5-)(<0)i=i  + 

*  oo  oo 

=  (£-)(«/>),'=.  +  (4.27) 

1  OO  1  oo 

Thb  is  ,  *  -  “pure”  extrapolation  on  the  density.  Using  (<5u),~o  —  [bT)j- o  -  0  in  equations 
(4.17)  and  (4.18)  leads  directly  to  a  —  1  and  (3  =  ^-1.  If  equation  (4.27)  were  a  pure 
extrapolation  on  p,  then  from  the  relation 

Icoobp  —  R\  +  R 2  +  2R3  (4.28) 

(see  equations  (4.9)  to  (4.11))  it  would  follow  that  o  =  1  and  t  =  2.  However,  it  should 
be  noted  that  while  6u  does  not  appear  in  equation  (4.27),  and  hence  o  —  1  in  (4.19), 
e  must  have  a  value  which  depends  on  the  solution  inside  the  domain.  Therefore,  it  is 
necessary  to  seek  solutions  to  the  nonlinear  system  (4.21)  and  (4.25)  for  a  given  A,  such 
that  |  z  |—  1,  z  7^  1,  |  ks  |<  1,  and  ao  =1.  The  solution  that  satisfies  z  —  e'e(6  ^  o), 

|  ks  |<  l(s  =  1,2,3)  with  ao  —  1  then  yields  the  value  of  c/3  —  (7  —  l)e.  The  numerical 
solution  of  equations  (4.21)  and  (4.25)  can  be  sensitive  near  |  ks  j=  1  so  double  precision 
was  used  in  the  numerical  solutions  of  these  equations.  For  a  wide  range  of  (and  hence 
A),  it  was  found  that  e/3  =  —1.24  ±  0.2,  giving  confidence  in  the  model2.  The  A  used  in 

2In  fact,  it  caji  be  easily  shown  that  the  model  boundary  conditions  (4.27)-(4.29),  with  a  =  o  =  1, 
(8  =  ^  —  1,  and  c  found  from  the  solution  to  (4.21)  to  (4.25),  correspond  to  a  primitive-variable  boundary 
condition  of  the  form  (p)J=o  =  (p)>=  +  a(^rm«i)2+™  where  a  =  0(1)  and  m  >  0.  Thus,  the  analysis 

is  performed  on  a  boundary  condition  wKich,  within  the  second-order  accuracy  of  the  finite-difference  code, 
cannot  be  numerically  distinguished  from  the  one  actually  imposed  i.e.,  equation  (4.26). 
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the  theoretical  prediction,  A  =  (u  +  c)^-,  was  based  on  the  actual  time  step  used  in  the 
Naviei-Stokes  codes  and  Ai  was  taken  as  the  radial  cell  size  (Armil/)  nearest  lo  the  inflow 
boundary. 


The  present  theoretical  model  most  closely  corresponds  to  the  Navier-Stokes  code  with 
the  finite-difference  method  in  the  radial  direction  and  the  Fourier  collocation  method 
in  the  circumferential  direction.  A  comparison  of  the  theoretical  Strouhal  number  from 
the  one-dimensional  model  with  the  one  obtained  from  this  Navier-Stokes  code  is  given 
in  Table  2.  In  view  of  the  fact  that  the  analysis  is  based  on  a  one-dimensional  model, 
the  agreement  between  the  values  of  St2  predicted  by  the  model  and  the  computed  values 
is  reasonably  good.  Furthermore,  the  model  predicts  the  trend  of  decreasing  St2  with 
increasing  domain  size. 


N 

A  T 

'  max 

Atcodc 

^ model 

(^2)  code 

— 

-3 

J 

100 

.41 

.0053 

.0452 

.0244 

.0296 

100 

.46 

.0053 

.0403 

.0225 

.0264 

Table  2:  Comparison  of  secondary  frequencies  predicted  by  the  analysis  and  computed 
using  the  mixed  spectral-finite-difference  code.  Mro  -0.4. 

Finally,  it  is  seen  that  if  the  characteristic  boundary  conditions,  (i.e.,  n  (1  0  in 

equations  (4.17)  to  (4.19))  are  used,  then  the  only  solution  to  (4.21)  through  (4.25)  is 
A 2  —  A3  --  0,  Ki  — -  1,  and  z  —  1.  It  can  be  shown  by  perturbation  analysis  that  it  is  not 
a  generalized  eigenvalue.  Thus,  any  perturbation  from  the  boundary,  R\  /li/Cj-s"  Aj 
will  remain  constant  and  small. 

5.  Concluding  remarks 

A  numerical  study  has  been  conducted  using  three  different  codes  to  compute  the 
unsteady  flow  about  a  circular  cylinder  in  low-speed  flow  at  low  Reynolds  numbers.  One 
of  these  codes  used  a  finite-difference  method  to  solve  the  two-dimensional  t  ime-dependent 
compressible  Navier-Stokes  equations,  the  second  used  spectral  techniques,  and  t  tie  third 
used  a  combination  of  these  two  methods.  With  stable  and  consistent  boundary  conditions, 
all  of  these  methods  were  able  to  predict  accurately  the  major  features  of  the  flow  such 
as  the  vortex-shedding  frequency.  However,  it  was  found  that  certain  far-field  boundary 
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conditions  which  used  extrapolation  of  the  primitive  variables  introduced  an  additional 
temporal  frequency  into  the  solution.  By  coincidence,  the  value  of  this  frequency  was 
very  close  to  that  found  experimentally  by  Sreenivasan  (1985).  The  use  of  characteristic 
variables  in  the  far-field  boundary  conditions  eliminated  this  frequency  from  the  solution. 
An  explanation  of  this  behavior  was  provided  using  an  analysis  based  on  a  simplified  model 
for  the  boundary  conditions.  These  results  illustrate  the  fact  that  great  care  must  be 
taken  in  interpreting  the  results  from  numerical  simulations,  and  that  while  the  secondary 
frequencies  are  spurious  in  the  infinite-domain  case,  the  possibility  exists  that  they  may 
be  found  in  a  confined  space  such  as  a  wind  tunnel.  In  fact,  Sreenivasan  3  found  that 
the  secondary  frequency  disappeared  when  the  tests  were  performed  in  a  wind  tunnel 
where  the  distance  from  the  cylinder  to  the  upper  and  lower  tunnel  walls  was  250  cylinder 
diameters.  This  distance  is  about  seven  times  that  used  in  the  original  tests.  Insight  into 
the  presence  of  the  multiple  frequencies  in  the  original  experimental  results  would  require 
three-dimensional  computations  which  included  the  wind-tunnel  walls. 
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PHYSICAL  SOLl/TION  DOMAIN 


FIGURE  1.  Solution  domain  In  physical  plane. 
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(a.)  Primitive-variable  boundary  conditions. 


(b.)  Characteristic-variable  boundary  conditions. 


FIGURE  2.  Lines  of  constant  pressure  (-M.  Finite-difference  code.  Contour  increment 

'  Poo 

=0.01.  My,  =  0.1,  =  80. 


FIGURE  3.  Linos  of  constant  pressure  (^).  Fully-spectral  code.  Contour  increment 
=0.01.  Moo  -  0.4,  Re =  80. 
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Primitive-variable  boundary  conditions. 


step  number 

Characteristic-variable  boundary  conditions. 


(a.)  Time  history  of  pressure. 


log(psd) 


Primitive-variable  boundary  conditions 
Characteristic-variable  boundary  conditions 
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(b.)  Frequency  spectrum  of  pressure. 


FIGURE  4.  Effect  of  boundary  conditions  on  pressure  in  wake  region  at  location  10 
cylinder  diameters  downstream  of  cylinder  and  one  diameter  above  wake  centerline.  Finite- 
difference  code,  psd  =  power  spectral  density. 
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(a.)  Time  history  of  pressure. 


-1 


(b.)  Frequency  spectrum  of  pressure. 


FIGURE  5.  Computer!  pressure  in  wake  region  at  location  10  cylinder  diameters  down¬ 
stream  of  cylinder  and  one  diameter  above  wake  centerline.  Fully-spectral  code,  psd  —  power 
spectral  density. 


24 


NASA 


Report  Documentation  Page 


2  Government  Accession  No 


1  Report  No 

NASA  CR- 181998 
] CASK  Report  No.  90-16 
4  T itle  and  Subtitle 


SECONDARY  FREQUENCIES  IN  THE  WAKE  OF  A  CIRCULAR 
CYLINDER  WITH  VORTEX  SHEDDING 


3.  Recipient's  Catalog  No. 


5  Report  Date 

February  1990 

6  Performing  Organization  Code 


Saul  S.  Abarbanel  James  C.  Townsend 

Wai  Sun  Don 
David  Gottlieb 
David  H.  Rudy 

Performing  Organization  Name  and  Address 

Institute  for  Computer  Applications  in  Science 
and  Engineering 

Mail  Stop  132C,  NASA  Langley  Research  Center 
Hampton,  V A  23665-5225 

Sponsoring  Agency  Name  and  Address 

National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Hampton,  VA  23665-5225 


8.  Performing  Organization  Report  No 

90-16 

10.  Work  Unit  No 

505-90-21-01 

11  Contract  ot  Grant  No 

NAS  1-1 8 107 
NASI -18*05 

13.  Type  of  Report  and  Period  Covered 
Contractor  Report 
14  Sponsoring  Agency  Code 


15  Supplementary  Notes 

Langley  Technical  Monitor: 
Richard  W.  Barnwell 


Submitted  to  Journal  of 
Fluid  Mechanics 


Final  Report 


16  Abstract 

A  detaiicd  numerical  study  of  two-dimensional  flow  past  a  circular  cylinder 
at  moderately  low  Reynolds  numbers  has  been  conducted  using  three  different  numer¬ 
ical  algorithms  for  solving  tlte  time-dependent  compressible  Navier-Stokes  equa¬ 
tions.  It  was  found  that  if  the  algorithm  and  associated  boundary  conditions  were 
consistent  and  stable,  then  the  major  features  of  the  unsteady  wake  were  well- 
predicted.  However,  it  was  also  found  that  even  stable  and  consistent  boundary 
conditions  could  introduce  additional  periodic  phenomena  reminiscent  of  the  type 
seen  in  previous  wind-tunnel  experiments.  However,  these  additional  frequencies 
were  eliminated  by  formulating  the  boundary  conditions  in  terms  of  the  character¬ 
istic  variables.  An  analysis  based  on  a  simplified  model  provides  an  explanation 
for  this  behavior. 


J  17  Key  Words  (Suggested  by  Authorlsll 

j  Navier-Stokes,  Low  Reynolds  numbers, 
j  cylinder,  CFD,  shedding  frequencies, 
boundary  conditions 


18  Distribution  Slafement 

02  -  Aerodynamics 
34  -  Fluid  Mechanics 


19  Security  Classif  (of  this  report) 
Unc lass  1 f i  ed 


_ _ Unclassified  -  Unlimi t ed 

[20  Security  Classif  (of  this  page)  21  No  of  pages 

!  Unclassified  26 


22  Price 


NASA  FORM  1626  OCT  86 


NASA  IahrIpv,  1990 


